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ABSTRACT 

In the paradigm of hierarchical galaxy formation, luminous radio galaxies mark mass assembly 
peaks that should contain clusters of galaxies. Observations of the z — 1.53 quasar 3C 270.1 with 
the Spitzer Space Telescope at 3.6-24 /jm and with the 6.5-m MMT in the z'- and Y"-bands allow 
detection of potential cluster members via photometric redshifts. Compared with nearby control 
fields, there is an excess of ~11 extremely red objects (EROs) at 1.33 < z p hot < 1.73, consistent 
with a proto-cluster around the quasar. The spectral energy distributions (SEDs) of 3/4 of the 
EROs are better fitted with passive elliptical galaxies than with dust-reddened starbursts, and of 
four sources well-detected on an archival HST snapshot image, all have undisturbed morphologies. 
However, one ERO, not covered by the HST image, is a double source with 0'.'8 separation on 
the z' image and a marginal (2er) 24 /zm detection indicating a dust-enshrouded starburst. The 
EROs are more luminous than L* (H — —23.6 AB mag at z ~ 1.5). 

Subject headings: Galaxies: clustering - quasars: general - Infrared: galaxies 



1. Introduction 

Groups and clusters of galaxies are key targets 
for studying the evolution of galaxies and large- 
scale structures. In the local universe thousands 
of clusters are known, but at high redshift (z > 1) 
only few clusters have been identified due to the 
challenge of discriminating a distant cluster from 
the abundant foreground galaxies. 

A few high-redshift clusters have been found via 
X-ray observations of a hot intra-cluster medium 
up to z — 1.45 (Mullis et al 2005, Bremer et al. 
2006, Stanford et al. 2006, Fang et al. 2007). Clus- 
ters up to this redshift range have also been found 
by searching for galaxy overdensities in wide-area 
optical-infrared surveys, for instance the Spitzer- 
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IRAC shallow survey (Stanford et al. 2005, Eisen- 
hardt et al. 2008), the UKIDSS ultra-deep survey 
(van Breukelen et al. 2006), HIROCS/COSMOS 
(Zatloukal et al. 2007), and the Spitzer First Look 
Survey (Muzzin et al. 2008). 

An alternative way to find high-redshift clusters 
is to use radio sources as signposts. Luminous ra- 
dio sources occur in massive galaxies (> 10 11 M Q , 
Seymour et al. 2007) , which according to the hier- 
archical formation model represent mass assembly 
peaks of the early universe. As such, they are 
expected to trace overdense regions and therefore 
cluster formation sites. The invaluable advantage 
of searching for clusters around radio sources is 
that one knows the redshift to search and there- 
fore one can limit the observational effort. 

Searches around very-high-redshift (2 < z < 5) 
radio sources have found an excess of Lyman a- 
cmitting galaxies (Kurk et al. 2000, Miley et al. 
2004, Venemans et al. 2007). These proto-clusters 
show a rather large projected extent and only a 
loose concentration, consistent with clusters that 
are forming but not yet virialized (Interna et al. 
2006, Zirm et al. 2008). 

Based on small (<2') maps, tentative (2-3cx) 



evidence for an overdensity of extremely red ob- 
jects (EROs, R — K > 5) around high-redshift 
(1 < z < 3) radio sources has been reported (Best 
2000, Sanchez & Gonzalez-Serrano 2002, Stern 
et al. 2003, Toft et al. 2003, Wold et al. 2003, 
Zirm et al. 2008). These EROs could be high- 
redshift passive ellipticals or dust-reddened star- 
forming galaxies. This ambiguity leaves open the 
question of whether the giant ellipticals were al- 
ready in place at these redshifts, implying a for- 
mation redshift larger still, or whether the EROs 
are starbursts and we are probing the epoch at 
which the cluster ellipticals formed. 

Using the Spitzer Space Telescope (Werner et al. 
2004), we are performing a 3.6-24 yum survey of 
the complete 3CR sample of radio galaxies and 
quasars at 1 < z < 2.5 (Haas et al. 2008). The 
large (>5') Spitzer arrays enable us to explore 
also the environment of the high-redshift radio 
sources. The 3.6-24 /im wavelength range mea- 
sures the redshifted stellar flux peak and the be- 
ginning of the Rayleigh-Jeans tail, while dusty 
starbursts may in principle show enhanced 24 /im 
emission. This makes the survey particularly well- 
suited to identify massive ellipticals and dusty 
starbursts. This paper presents first results from 
a pilot study of the quasar 3C 270.1 at z = 1.53, 
combining our Spitzer data with z' - and V-band 
maps. Distances and ages in this paper are based 
on a ACDM cosmology with Hq = 71 km/s/Mpc, 
Q m = 0.27 and 0\ = 0.73. 

2. Observations and data 

The Spitzer maps of 3C 270.1 were obtained us- 
ing the instruments IRAC (3.6-8.0 /im, Fazio et al. 
2004) and MIPS (24 /im, Rieke et al. 2004). On- 
source exposure times were 4x30s (each IRAC 
band) and 10 x 10s (MIPS). The dithered maps 
(PID 40072, PI G. Fazio) cover about 4' FOV with 
full depth. In addition, comparison fields offset 6.'5 
from the main field were automatically obtained in 
the IRAC bands, at 3.6 /im and 5.8 /im northeast 
and at 4.5 /im and 8.0 /im southwest of the central 
field containing the quasar. The maps allow us 
to detect sources down to ^4uJy (3a) at 3.6 and 
4.5 fim and 100 /LtJy at 24/im0 



1 We have also obtained a 16 fim map with the IRS peak-up 
array (4 x14 s, Houck et al. 2004) which has been useful for 
the quasar SED (Haas et al. 2008). Because of the lim- 



Complementary z'- and V-band images were 
obtained at the 6.5-m MMT using MegaCam (30' 
FOV, McLeod et al. 2006) and SWIRC (5' FOV, 
Brown et al. 2008). The z'- and V-bands (cen- 
tered at 0.89 and 1.05 /im, respectively) bracket 
the 4000 A break at redshift z = 1.53. Sensitiv- 
ities better than 1 /iJy were achieved at z'- and 
V-band under good seeing conditions (FWHM < 
1") with exposure times of 40 and 90 min, re- 
spectively. The z'-band image covers the IRAC 
comparison fields, but the Y image does not. At 
Galactic latitude b = 80?6, foreground extinction 
is negligible (Eb-v — 0.012 mag). 

The images were reduced using interactive anal- 
ysis tools. The IRAC mosaics were corrected for 
the residual images that arise from prior observa- 
tions of bright sources by making object-masked, 
median-stacked coadds with all the science frames 
and then subtracting them from the science frames 
prior to final coadding. Thus, we can be sure 
there are no faint spurious signals in the IRAC mo- 
saics arising from prior observations. For IRAC, 
we used the basic calibrated data products (BCD, 
version S16) and coadded them to 0'.'869 pixels 
using version 4.1.2 of IRACProc (Schuster et al. 
2006). This optimally handles the slightly under- 
sampled IRAC PSF in order to assure accurate 
point-source photometry. For MIPS, we used cus- 
tom routines to modify the version S17 BCD files 
to remove instrumental artifacts (e.g., residual im- 
ages) before shifting and coadding to create the 
final mosaics. The MegaCam and SWIRC frames 
were reduced using standard procedure^. 

The 3.6 /xm image shows more sources than any 
other, and we used it as the source detection im- 
age for constructing a multi-band catalog. In the 
central 4' x 4' with optimal depth, 274 objects 
are detected at 3.6 /im, with 272/202/212 of them 
being also detected at 4.5/5.8/8.0/im At shorter 
wavelength, 184 (150) of the 3.6 /im sources are 
detected at z' (Y), and 143 at both z' and Y . 
The sources were extracted and matched using 
the SEXTRACTOR tool (Bertin & Arnouts 1996). 
Details are described in Appendix IA~1 The photo- 
metric uncertainties are typically less than 10% 



ited sensitivity (100 /xjy) and small size (l.'5x2'), this map 
does not provide useful constraints for the cluster search. 
Therefore, the 16 fim data are not further considered here. 

2 See http:/ /www. cfa.harvard.edu/ ~mashby/megacam/megacam_frames.html 
for MegaCam. 
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but increase for the faintest sources. 

The region of 3C 270.1 has also been mapped 
with the Sloan Digital Sky Survey (SDSS, DR6). 
About 35 relatively bright SDSS sources are iden- 
tified with foreground galaxies and stars. These 
objects are excluded from the following discussion. 
In addition, we found an archival Hubble Space 
Telescope WFPC2 snapshot image of 3C 270.1 in 
the F702W passband (600 s exposure time, PI W. 
Sparks, Lehnert et al. 1999). Because of the lim- 
ited spatial coverage and sensitivity (0.6 /iJy), we 
did not use the F702W photometry for finding po- 
tential cluster members. However, the HST image 
enables us to inspect the morphology of a few of 
the sources to determine whether they are point- 
like or extended and to look for rest-frame UV 
excesses. 

3. Results and discussion 
3.1. Cluster galaxy candidates 

Cluster galaxies around 3C 270.1 should lie at 
the redshift of the quasar. We determined photo- 
metric redshifts z p hot by fitting the observed spec- 
tral energy distributions (SEDs) of the 184 galax- 
ies detected at z' with two basic templates, an 
elliptical galaxy and a dusty starburst galaxyH 
For the elliptical galaxy, we used NGC 221, which 
has a strong 4000 A break. The NGC 221 spec- 
tral template from Kinney et al. (1996), covering 
the rest-frame wavelength range 3000-8000 A, was 
smoothly extrapolated to longer wavelengths by a 
4000 K blackbody. For the starburst galaxy tem- 
plate, we used the ultra-luminous infrared galaxy 
Arp220 with photometry from the NASA Ex- 
tragalactic Database (NED) and the Sloan Dig- 
ital Sky Survey (SDSS, DR6) and mid-IR spec- 
tra from Spitzer/IRS. We also tested the dust- 
reddened star- forming galaxy template M82, but 
the results were similar to those for Arp 220 (dif- 
ferences in Zphot less than Az = 0.1), and therefore 
we present only the Arp 220 results. 

The accuracy of photometric redshifts can be 
estimated, for instance, from the Spitzer Wide- 
Area Infrared Extragalactic Survey (SWIRE), 
which has spectroscopic redshifts available for 



3 We also tried type-l/-2 AGN templates, but none of the 
SEDs is consistent with such templates at redshift z ~ 1.5 
except 3C 270.1 itself. 



many of its sources. In SWIRE, 7 filters (5 op- 
tical and 2 IRAC) were sufficient to discriminate 
between 8 templates (ranging from blue to red 
galaxy types) and to determine z p h ot with an rms 
of Az/(1 + z) = 3.5% (Rowan-Robinson et al. 
2008). We have fewer filters than SWIRE, but 
our task is easier because we only have to deter- 
mine whether the photometry is consistent with 
z = 1.53 or not. Furthermore, we consider only 
extremely red sources and only two templates and 
therefore suggest that we can reach accuracy sim- 
ilar to that of SWIRE. At z — 1.5, an accuracy 
(rms) of about 4% corresponds to Az = 0.1. This 
is much larger than the expected redshift disper- 
sion within a cluster (Az < 0.01). While the SEDs 
can provide cluster galaxy candidates within ap- 
propriate redshift bins around zqsoi confirmation 
of a cluster will require spectroscopic redshifts of 
at least a sample of the candidates. 

Figures [1] and [5] show examples of SEDs and 
template fits to cluster galaxy candidates. The 
most striking SED features are the steep rise from 
z'- to y-band and beyond and the decline between 
3.6 and 8.0 fim. The slopes of these features de- 
termine the redshift. We performed the template 
fits over a grid in redshift (dz = 0.01) and inten- 
sity. For a red SED with F(3.6)/F(z') w 30, the 
chi square contour plots exhibit a sharp minimum, 
suggesting that for a given template the achievable 
accuracy of z p hot is about Az = 0.1. As illus- 
trated in Fig. [21 even for SEDs with fewer than 
six data points the accuracy of z p hot should be 
not worse than Az < 0.20. Therefore, we have 
chosen as cluster galaxy candidates the 29 objects 
for which either the elliptical or the starburst tem- 
plate yields 1.33 < z p hot < 1-73. The basic con- 
clusions on the clustering of EROs around 3C270.1 
remain unchanged for other ranges |Az| between 
0.15 and 0.25. 

Usually the 0.9-8.0 /im SEDs can be fitted very 
well with both the elliptical and the starburst tem- 
plate. Both templates are extremely red, and 
therefore no other galaxy type is likely to fit the 
SEDs. The examples in Figures Q] and [2] illus- 
trate the ability to identify the two galaxy types 
at z ~ 1.53. The resulting fit parameters for the 
29 cluster candidates are listed in Table [T] The 
photometric redshifts differ systematically for the 
elliptical and starburst templates. On average, 
z e ii is higher by 0.27 ± 0.1 than zsb- Among 
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Fig. 1. — Observed spectral energy distributions 
for two good examples of the 29 cluster candi- 
dates. The IRAC and MMT data are marked with 
filled circles and la error bars. The MIPS 24 /jm 
data point is marked with a filled circle, too; in 
the upper panel (elliptical source) it is a 3cr upper 
limit, in the lower panel (dusty starburst source) 
it is a 2cr detection also visible on the map. The 
horizontal bar indicates the 24 /im pass band for 
comparison with the silicate absorption feature. 
HST photometry is marked with an open circle; 
it is not available for the dusty starburst source. 
Solid lines show the elliptical galaxy (NGC 221) 
fit for the source in the upper panel and the dusty 
starburst (Arp 220) fit for the galaxy in the lower 
panel. Dotted lines show the alternative template 
for each galaxy. The upper and lower panel shows 
object numbers 23 and 3, respectively, as listed in 
Tabled]). 

the 29 cluster candidates, the elliptical template 
yields better fits (smaller chi square values) for 
22 sources (75%, Table [TJ. All of these sources 
lie in the redshift bin 1.33 < zsb < 1.73. Among 
the remaining 7 cluster candidates, those 3 sources 
which favor the starburst fit lie in the redshift 
bin 1.33 < zsb < 1-73, too. These sources have 
Xsb < 0-3 ' XeU (object numbers 3, 5 and 29 in Ta- 
ble [T]). The remaining 4 sources (object numbers 
11, 14, 17 and 28) are not detected at 5.8 and 8.0 
fim, and they do not clearly favor the starburst 
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Fig. 2. — Observed spectral energy distributions 
for two cluster candidates with poorer data qual- 
ity. The IRAC and MMT data are marked with 
filled circles and la error bars. The sources are de- 
tected only in the z'-band (in the lower panel also 
in the V-band) and at 3.6 and 4.5 /im, but the up- 
per limits at 5.8 and 8.0 /im help to constrain the 
redshift fits. As in Fig. [TJ the upper panel shows 
a source (object 11) which is preferably fit by the 
elliptical template and the lower panel one (29) fit 
by the starburst template. The dotted and dashed 
lines show the preferred templates at different red- 
shifts (indicated in the figure), showing that the 
accuracy of the photometric redshifts should be 
dz < 0.2. 

template, so that - as ellipticals - they are cluster 
galaxy candidates. 

The cluster galaxy candidates occupy a limited 
range of z' — [3.6] color as shown in Figure [31 
While the candidates were determined from multi- 
band SED fitting (usually 4-6 bands), selecting a 
color range of 15 < F(3.6)/F(z') < 52 would have 
given nearly the same sample. The sources red- 
der than this color range are probably at larger 
redshift (z > 1.8), while bluer sources are consis- 
tent with being either foreground objects or unred- 
dened star-forming galaxies at the quasar redshift. 
Our EROs with F(3.6)/F(z') > 15 (correspond- 
ing to z' — [3.6] > 5 mag in the Vega system) are 
likely to also obey the standard ERO definition of 
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R- K > 5 mag (cf. Wilson et al. 2004). 
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Fig. 3. — Color-magnitude diagram F(3.6)/ F(z') 
versus -F(3.6), for the quasar field (upper panel) 
and control field 1 (lower panel). Sources with 
Zphot = 1-53 ± 0.20 (SED determined candidates 
marked with circles) concentrate in a distinct color 
range 15 < F(3.G)/F(z') < 52 indicated by the 
horizontal dotted lines. The number of color de- 
termined candidates is also given. 



3.2. Sky distribution and surface density 

Figure [4] shows the sky distribution of the clus- 
ter candidates. Spreading over a diameter of more 
than 2', they form at most a loose concentra- 
tion around the quasar, indicating a proto-cluster 
rather than a virialised, concentrated system. The 
peak density appears to lie ~30" to the east of 
3C 270.1, but it is difficult to be precise with so 
few galaxies. The 3C 270.1 proto-cluster is larger 
in angular size than the X-ray clusters at z = 1.45 
and z = 1.22 found by Mullis et al. (2005) and Bre- 
mer et al. (2006). These clusters have an extent of 
less than 1' in both their galaxy distributions and 
their X-ray sizes. 
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Fig. 4. — Sky distribution of the 29 candidate clus- 
ter galaxies around the quasar 3C 270.1 (marked 
with a cross). The triangle shows the apparent 
centroid of the cluster galaxy distribution. The 
star in the southwest marks the starburst candi- 
date whose SED is shown in Fig.[T] The solid lines 
surround the areas covered by IRAC and HST 
frames. The dotted circles of radius 50", 100", 
and 150" outline the areas considered in Fig.EJ 
they are centred around the quasar. At z = 1.53, 
50" corresponds to 427 comoving kpc. 

We have also analyzed the two comparison 
fields taken with IRAC and covered by the Mega- 
Cam z' image (but not the SWIRC Y image). De- 
tails of the analysis are described in Appendix [B] 
The surface density of possible z = 1.53 ± 0.20 
galaxies is less than about 18 objects in each con- 
trol field of area 15.2 arcmin 2 (1.2 per arcmin 2 ). 
Thus, in the redshift range z — 1.53 ± 0.20, the 
central field surrounding the quasar shows an ex- 
cess of at least about 29 — 18 — 11 sources, i.e., 
60% over the comparison fields. 

Figure [5] shows the radial surface density plot of 
the 29 EROs in the redshift range z = 1.53 ±0.20. 
The surface density peaks inside the central 50" 
radius and declines steadily with increasing dis- 
tance down to the surface density of the control 
fields. This provides further evidence that there is 
an excess of z p hot ~ 1-53 galaxies near 3C 270.1. 
The radial overdensity is also present when using 
circles around the centroid east of the quasar po- 
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Fig. 5. — Surface density of the 29 cluster galaxy 
candidates versus projected distance from the 
quasar 3C270.1 (solid line through fat dots with 
Poisson error bars). The radial bins centred 
around the quasar are outlined in Fig. [H and the 
surface density of the outermost annulus has been 
corrected for the area not covered by IRAC. For 
comparison, the long-dashed histogram shows the 
surface density of cluster galaxy candidates versus 
projected distance from centroid (the triangle in 
Fig. 0| . The dotted line indicates the mean sur- 
face density in the two control fields. 



sition. 

Existing X-ray data on 3C 270.1 are inconclu- 
sive about the presence of a cluster. The quasar 
itself is clearly detected as a point source in our 
10 ks Chandra data (Wilkes et a!., in prepara- 
tion). Weak extended X-ray emission (<20") is 
also present, but most of the counts come from 
the position of the southern radio lobe or from 
between the northern radio lobe and the quasar. 
ROSAT data (19.3 ks, 1993 May) show a strong 
detection of 3C 270.1, but at ROSAT resolution 
not only is there no way to separate cluster gas 
emission from quasar emission, but the emission is 
also heavily blended with that of unrelated QSO 
B1223+338B located 0f8 awayH We found no 
XMM observations. Longer Chandra exposures 



4 Because of its lower rcdshift, z = 1.038, this QSO does not 
affect the clustering evidence discussed here. 



will be needed for a definitive detection of any X- 
ray emission from cluster gas. 

3.3. Nature of cluster galaxy candidates 

The cluster galaxy candidates have an abso- 
lute rest-frame magnitude H « K < —23.8 (AB 
system). An L* galaxy at z ~ 1.5 has H w 
K = —23.6 (as determined from observations at 
1 < z < 1.3 by De Propris et al. 2007). Thus 
only the most luminous cluster galaxy candidates 
are detected on our maps. The galaxies are either 
giant ellipticals or dust-reddened starbursts. If gi- 
ant ellipticals, most of their stellar mass has been 
formed at even higher redshift, but they may har- 
bor some ongoing star-formation — it would be 
relatively weak with respect to the already exist- 
ing stellar mass but could show up at rest UV 
wavelengths. If dusty starbursts, they would have 
to be at least three times more luminous in the 
near-infrared than the ULIRG Arp220, indicating 
that a large stellar mass has already formed during 
an earlier episode. 

The distinction between ellipticals and dusty 
starbursts is difficult to make at rest wavelengths 
shorter than 1 /im using photometry only, al- 
beit possible with adequate wavelength coverage 
(Pozzetti & Mannucci 2000). At z = 1.5, the six 
data points of our SEDs cover about rest wave- 
length 0.3-3 /im. As mentioned in Sect. 13. 1\ 75% 
of the SEDs are better fitted by the elliptical than 
the starburst template. In principle, additional 
sensitive photometry at rest wavelengths longer 
than 10 /im could improve the starburst-elliptical 
distinction (Stern et al. 2006). At z = 1.53, unfor- 
tunately, the 9.7/tm silicate absorption enters the 
24 /^m band, reducing the potential to detect the 
powerful MIR emission of starbursts. Apart from 
3C 270.1 itself, there is only one cluster candidate 
with a marginal (2a) detection at 24 /tm (Fig. [2]). 
Its SED provides evidence for a dust-enshrouded 
starburst. While only elongated on the F-band 
image, on the z'-band image this galaxy appears 
as a double source with 0'.'8 (^7kpc) separation; 
for determining z p hot we have used the combined 
z'-band photometry. This source is not covered 
by the HST image (Fig. g]). 

Because the z'-band image is limited in spatial 
resolution, we have also inspected the morphology 
of the sources seen on the HST snapshot image. 
The HST frame covers 15 potential cluster galax- 
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ies (Fig. [4|, 9 of which are detected. All 9 are 
extended, ruling out the (remote) possibility that 
they might have been brown dwarfs. Figure [6] 
shows the four sources detected in the F702W 
filter with photometric accuracy better than 5<r. 
The SEDs of these galaxies are best fit with the 
elliptical template. They have a regular shape, 
and none of them shows a peculiar morphology. 
This argues against starburst galaxy pairs and in 
favor of an evolved elliptical population. The re- 
maining five of the 9 HST-detected sources are too 
faint to draw stringent morphological conclusions. 
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Fig. 6.— HST F702W images of four cluster 
galaxy candidates (object numbers 9, 12, 14 and 
15 in Table [T]). The contours are linearly spaced 
in steps of 10% of the peak flux value (0%, 10%, 
... 90%). Arrows indicate the orientation of each 
panel. 

The HST F702W photometry is of limited use 
because of the short exposure time, but six of the 
nine detected sources appear to show a rest-frame 
UV excess above the elliptical template (Fig. [2]). 
This indicates some level of ongoing star forma- 
tion activity. However, at rest wavelengths around 
1 /zm this activity appears to be outshone by large 
numbers of already-evolved stars. 

Having so many as ~11 galaxies brighter than 
L* would be very high for a local cluster, but pas- 
sive evolution will cause these galaxies to become 



fainter by something like 3 ma^ by z — to 
become comparable to today's cluster ellipticals. 
Thus the photometric data are compatible with 
the existence of a cluster. 

4. Conclusions 

High redshift radio sources are predicted to 
serve as signposts for galaxy clusters in the early 
universe. In a pilot study using sensitive, large 
maps from the Spitzer Space Telescope and the 
6.5-m MMT, we find a clustering of red galaxies 
around the z — 1.53 quasar 3C 270.1: 

1) Photometric redshifts identified 29 galaxies 
consistent with the quasar's redshift. While 
moderate clustering is evident in both the 
radial surface density of these objects and in 
their numbers with respect to nearby control 
fields, the sky distribution of the EROs sug- 
gests a loose, not-yet-virialized protocluster. 
Deep, high- angular-resolution X-ray obser- 
vations are required to see whether hot intra- 
cluster gas is present or not. 

2) The cluster-candidate galaxies have z' — [3.6] 
colors indicative of extremely red objects 
(EROs). They may be either dusty star- 
bursts or elliptical galaxies with little ongo- 
ing star formation. The SED fits favor the 
majority of the potential cluster members 
being passive ellipticals, and this is consis- 
tent with the four objects with good HST 
images having undisturbed morphology. In 
some cases the HST photometry reveals a 
rest-frame UV excess over the elliptical tem- 
plate, consistent with modest on-going star 
formation and/or nuclear activity in these 
galaxies, but they must have formed the bulk 
of their stars at redshifts >1.53. 

This pilot study demonstrates that the Spitzer j 
IRAC maps provide an efficient way to search for 
clustering of red galaxies around high redshift ra- 
dio sources, but accurate redshifts and the nature 
of the galaxies have to be confirmed with addi- 
tional spectroscopy and/ or deep far-infrared imag- 
ing with the Herschel Space Observatory. The 



5 Passive evolution was estimated via a Starburst99 model 
(Vazquez & Leitherer 2005) with an instantaneous burst 
age 1 Gyr before z = 1.53 and thus an age of 10 Gyr at 
z = 0. 
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ongoing investigation of all 64 high-redshift 3CR 
sources will result in a homogeneous database of 
considerable cosmological impact. 
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Fig. 7. — Color-magnitude diagram F(4.5)/F(z') 
versus F(4.5) of the quasar field and control field 
2. top: quasar field with all sources detected at 

4.5 ^m using SExtractor in double-image mode at 

3.6 & 4.5 yum. middle: quasar field restricted to 
those sources detected at 4.5 ^im using SEctractor 
in single-image mode, as was done for the con- 
trol field 2. bottom: control field 2. Sources with 
Zphot = 1-53 ± 0.20 (SED determined candidates 
marked with circles) concentrate in a distinct color 
range 13 < F(4.5)/F(z') < 54 indicated by the 
horizontal dotted lines. The number of color de- 
termined candidates is also given. 
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Table 1: Fit parameters for the 29 cluster galaxy candidates in the quasar field. 



Object RA J2000 Dec J2000 n a z ell b X jn z SB c 



1 


12 


20 


27 


.90 


33 


43 


11 


.5 


3 


1 


.47 


0. 


.43 


1 


.16 


0.53 


2 


12 


20 


28 


.72 


33 


41 


51 


.9 


6 


1 


.50 


0. 


.06 


1 


.20 


0.22 


3 


12 


20 


29 


.20 


33 


42 


10 


.5 


6 


1 


.55 


0. 


.51 


1 


.35 


0.15 


4 


12 


20 


30 


.06 


33 


43 


18 


.1 


6 


1 


.47 


0. 


,19 


1 


.00 


0.85 


5 


12 


20 


30 


.66 


33 


44 


50 


.0 


4 


1 


.90 





.17 


1 


.60 


0.03 


6 


12 


20 


30 


.73 


33 


44 


39 


.4 


6 


1 


.53 


0. 


.12 


1 


.24 


0.35 


7 


12 


20 


31 


.18 


33 


43 


14 


.1 


6 


1 


.63 


0. 


.04 


1 


.47 


0.31 


8 


12 


20 


32 


.49 


33 


43 


20 


.3 


5 


1 


.50 


0, 


,29 


1 


.34 


0.39 


9 


12 


20 


33 


.50 


33 


42 


25 


.4 


5 


1 


.35 





,10 


1 


.11 


0.38 


10 


12 


20 


33 


.57 


33 


44 


08 


.3 


4 


1 


.42 


0. 


.13 


1 


.16 


0.18 


11 


12 


20 


33 


O CT 

.85 


33 


45 


37 


,5 


3 


1 


.53 


0. 


.07 


1 


riri 
,11 


0.03 


12 


12 


20 


33 


.90 


33 


42 


34 


.6 


4 


1 


.36 





,08 


1 


.17 


0.30 


13 


12 


20 


34 


.79 


33 


41 


55 


.4 


5 


1 


.47 


0. 


.10 


1 


.07 


0.27 


14 


12 


20 


34 


.95 


33 


42 


41 


.0 


4 


1 


.45 


0. 


,01 


1 


.20 


0.01 


15 


12 


20 


35 


.20 


33 


42 


09 


.6 


5 


1 


.35 





,12 





.96 


0.43 


16 


12 


20 


35 


.88 


33 


43 


32 


.2 


4 


1 


.59 





,07 


1 


.28 


0.10 


17 


12 


20 


36 


.50 


33 


43 


30 


.6 


4 


1 


.40 


0, 


.19 


1 


.20 


0.10 


18 


12 


20 


37 


.70 


33 


41 


18 


.2 


5 


1 


.54 


0, 


.34 


1 


.33 


0.37 


19 


12 


20 


38 


.52 


33 


43 


04 


.2 


4 


1 


.53 


0. 


.05 


1 


.25 


0.08 


20 


12 


20 


38 


.78 


33 


43 


50 


.2 


6 


1 


.68 


0. 


,32 


1 


.38 


0.53 


21 


12 


20 


39 


.16 


33 


44 


31 


.9 


4 


1 


.50 


0. 


,01 


1 


.27 


0.03 


22 


12 


20 


40 


.10 


33 


43 


20 


.5 


6 


1 


.50 


0. 


,23 


1 


.36 


0.35 


23 


12 


20 


40 


.49 


33 


43 


21 


.9 


6 


1 


.56 


0. 


.13 


1 


.04 


0.98 


24 


12 


20 


40 


.55 


33 


43 


17 


.3 


4 


1 


.63 


0. 


.03 


1 


.55 


0.06 


25 


12 


20 


40 


.62 


33 


43 


00 


.2 


4 


1 


.48 


0. 


,01 


1 


.25 


0.02 


26 


12 


20 


43 


.31 


33 


42 


51 


.2 


5 


1 


.50 





.44 


1 


.32 


0.66 


27 


12 


20 


43 


.32 


33 


44 


19 


.8 


5 


1 


.69 


0. 


.01 


1 


.37 


0.04 


28 


12 


20 


44 


.17 


33 


44 


03 


.9 


4 


1 


.56 


0. 


.11 


1 


.27 


0.07 


29 


12 


20 


45 


.13 


33 


43 


18 


.9 


4 


1 


.84 


0. 


.26 


1 


.53 


0.04 



"Number of data points (detections) used to calculate x 2 of the fit. 
kftedshift of the fit using the elliptical template NGC 221 
c Rcdshift of the fit using the starburst template Arp 220 
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A. Source extraction and photometry 

For source extraction and photometry, we used the SExtractor tooH (Bertin & Arnouts 1996, Version 
2.5). This tool is able to identify sources via a tree structure algorithm, provide their total magnitudes, 
and, if a saddle point is present in their profile, to deblend components via a multi-isophotal analysis. 
Photometry is derived using automatically adjusted apertures (option MAG-AUTO) following Kron's first 
moment algorithm (Kron 1980, Infante 1987). In case close neighbours are found, the components are 
iteratively deblended, with the photometric uncertainty of the faint companion reaching 25% in the worst 
cases (Beard et al. 1991). 

In a first step we adjusted the SExtractor parameters for each image (single-image mode). The parameters 
are, for instance, noise threshold and minimum number of connected pixels defining a source. The MMT z'- 
and Y-band images had a seeing of about 0'.'9 (FWHM), the IRAC PSF has FWHM between l'/4 at 3.6 /im 
and l'/7 at 8.0 fim. Typical Kron radii of the sources (containing 90% of the flux) are about 3-4 arcsec in 
all filters (except 24 /im), so that the PSF differences play a minor role (< 5%), and the z'- and y-band 
photometry — as computed by SExtractor — needs no aperture correction^ For the IRAC photometry, 
however, the absolute photometric calibration refers to an aperture of 12", so we applied standard aperture 
corrections (5— 30% jH to derive the final photometry of the sources, all of which are essentially unresolved. 
The very few sources detected at 24 /jm are well separated, so that clean photometry could be derived 
without need for aperture correction. 

In a second step, we ran SExtractor in double-image mode on the 4.5, 5.8, and 8.0 /im IRAC maps, using 
as primary the 3.6 /im map because it contains more sources than any other IRAC map. The double-image 
mode mode is able to identify, for instance, for a 3.6 /im detection also the corresponding 4.5 /im source, 
even if it was not detected in single-image mode. The double-image mode enabled us to detect essentially 
all 3.6 /im sources at 4.5 /im (272 of 274 sources). The detection rate at 5.8 and 8.0 /im was lower (202 and 
212 sources, respectively). We determined the 5.8 and 8.0 /im upper flux limits (~7/iJy) from the break of 
the flux histograms of the detections. In case of non-detections and clean regions we used this flux value as 
the 3(7 upper limit. 

SExtractor was run separately on the z 1 image, and the z' and IRAC source lists were matched with 
a search radius of 1.5". This should produce good identifications because no two IRAC sources are closer 
together than 3". To check further how far confusion plays a role, we counted the number of MMT sources 
within a radius of 3" around each IRAC source. By far most sources (>90%) are unambiguously matched to 
a single z' source. Of 274 IRAC sources in the quasar field, 18 (6.6%) have two z 1 counterpart, and one even 
has three. The control fields show about the same statistics: 9 double z' counterparts (4.4%) out of 205 IRAC 
sources in the 3.6 /im field and 10 double counterparts (6.0%) out of 167 IRAC sources in the 4.5 /tm field. 
Thus, the fraction of confused sources is small and comparable for control and quasar fields, and therefore 
we do not expect that the clustering results are significantly affected by confusion. To avoid any possible 
complication, we have omitted from the analysis all sources (4-7%) with two or more z' counterparts found 
by SExtractor. The starburst cluster candidate galaxy (object number 3 in Table [T]) remained in the sample 
because SExtractor treated it as a single source. Most likely the saddle between the two components is not 
deep enough for SExtractor to split them despite the source's visual appearance on the z' image as a 0"8 
double. 



6 http:/ /terapix. iap.fr/IMG/pdf/sextractor.pdf 

7 Our MMT z'-band photometry was calibrated with the z-band image of the Sloan Digital Sky Survey (SDSS, DR6), including a 
(very small) color term to account for the similar but nonetheless different bandpasses. Our Y-band photometry was calibrated 
with interpolated total magnitudes from SDSS and 2MASS. 

8 http: / /ssc. spitzer.caltech.edu/irac/dh/ 
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B. Comparison with the control fields 

The quasar field has photometry in six filters, two MMT and four IRAC bands. But the two control fields 
have photometry in only three filters, in z' and two IRAC bands: 3.6 and 5.8 /im in control field 1, 4.5 and 
8.0 jum in control field 2. This makes the uncertainties on the photometric redshifts larger, and therefore 
comparison within the same redshift bin appears not to be straightforward. Furthermore, in the central 
field, the source detection rate at 4.5 /zm is lower than at 3.6 /im, so that for a fair comparison the number 
of candidates in control field 2 has to be corrected for. In order to facilitate the comparison, we used the 
colors z 1 — [3.6] and z 1 — [4.5] and counted the number of candidates in the respective color bins. 

Control field 1: The colors of the elliptical and starburst template, when redshifted to z = 1.53, are 
F(3.6)/F(z') = 25 and 39, respectively. In the color range 15 < F(3.6)/F(z') < 52 we found 28 and 
17 sources for the quasar and the control field, respectively (Fig. [3|). With regard to the low number 
statistics, these values are consistent with the 29 and 16 sources in the quasar and the control field, having 
1.33 < z p hot < 1-73 (denoted as SED candidates in Fig. [3|). This leaves an excess of 11 colour- and 13 
SED-selected sources in the quasar field. 

Control field 2: The colors of the elliptical and starburst template, redshifted to z = 1.53, are F(4.5)/ F(z') 
= 24 and 43, respectively. Guided by these "expectation values" and the actual z' — [4.5] distribution of the 
candidates in the quasar field, we selected a color range of 13 < F(4.5)/ F(z') < 54 (Fig. top). Now we 
considered the fields, where the 4.5 /zm source list was created using the single-mode SExtractor option, i.e., 
not making use of the 3.6 /im information in the quasar field. We found 24 and 14 color-selected candidates 
for the quasar and the control field, respectively (Fig. middle and bottom). These values are consistent 
with the corresponding numbers for SED-selected candidates of 23 and 12 in the quasar and the control 
field, respectively. Because the quasar field contains more cluster galaxy candidates, 29 via color and 29 
via z p hot (Fig. [3 top), we have scaled up the number of candidates in the control field by the factors 29/24 
for the color-candidates and 29/23 for the SED-candidates. This results in 17 colour- and 15 SED-selected 
candidates. This leaves an excess of 12 (= 29 — 17) colour- and 14 (— 29 — 15) SED-selected sources in the 
quasar field. 

Combining the control field counts we estimate an excess of 11-12 colour- and 13-14 SED-selected candi- 
date cluster galaxies in the quasar field. 
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